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After inserting the heavy quark mass dependence into QCD partonic evolution equations, we 
determine the mean charged hadron multiplicity of jets produced in high energy collisions. We 
thereby extend the so-called dead cone effect to the phenomenology of multip article production in 
QCD jets and find that the average multiplicity of heavy-quark initiated jets decreases significantly 
as compared to the massless case, even taking into account the weak decay products of the leading 
primary quark. We emphasize the relevance of our study as a complementary check of ^-tagging 
techniques at hadron colliders like the Tevatron and the LHC. 
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1. Introduction 

High-p ( jets can be initiated either in a short-distance interaction among partons in high en- 
ergy collisions such as pp, pp, in the DIS e ± p, the e + e~ annihilation and via electroweak (or 
new physics) processes. One well-known example is given by the decay chain of the top quark 
t — > H + b, where the b quark should start a jet. Thus the ability to identify jets from the frag- 
mentation and hadronization of b quarks becomes very important for such Higgs boson searches. 
Needless to say, the relevance of ^-tagging extends over many other channels in the quest for new 
physics at hadron colliders. The experimental identification of ft-jets relies upon several of their 
properties in order to reject background, e.g. jets initiated by lighter quarks or gluons. First, the 
fragmentation is hard and the leading ^-hadron retains a large part of the original b quark momen- 
tum. In addition, the weak decay products may have a large transverse momentum with respect 
to the jet axis therefore allowing separation from the rest of the cascade particles. Lastly, the rel- 
atively long lifetime of Z?-hadrons leading to displaced vertices which can be identified by using 
well-known impact parameter techniques [1]. Still, a fraction of light jets could be mis-identified 
as ft-jets, especially at large transverse momentum of the jet. Now, let us point out that an essential 
difference between heavy and light quark jets results from kinematics constraints: the gluon radia- 
tion off a quark of mass m and energy E ^> m is suppressed inside a forward cone with an opening 
angle & m = m/E, the so-called dead cone phenomenon [2,3]. 

In this work, we compute the average (charged) multiplicity of a jet initiated by a heavy 
quark. For this purpose, we extend the modified leading logarithmic approximation (MLLA) 
evolution equations [4] to the case where the jet is initiated by a heavy (charm, bottom) quark. 
The average multiplicity of light quarks jets produced in high energy collisions can be written 
as N(Y) oc exp < J y(y)dy > (Y and y are defined in section ^, where y ~ /o + Ay is the anoma- 
lous dimension that accounts for soft and collinear gluons in the double logarithmic approximation 
(DLA) /o ~ y/oT s , in addition to hard collinear gluons Ay ~ a s or single logarithms (SLs), which 
better account for energy conservation and the running of the coupling constant a s [4]. 

In the present work we evaluate the mean multiplicity as a function of the mass of the heavy 
quark. Actually, the mass depends on the scale Q of the hard process in which the heavy quark 
participates. In the MS renormalization scheme, for example, the running mass becomes a func- 
tion of the strong coupling constant CC S (Q), and the pole mass defined through the renormalized 
heavy quark propagator (for a review see [5] and references therein). In order to assess the scheme- 
dependence of our calculations beyond leading-order we have considered a broad range of possible 
values for the charm and bottom masses: m c = 1 — 1.5 GeV and mi, = 3 — 5 GeV, respectively [5], 
Lastly, we will see that under the assumption of local parton hadron duality (LPHD) as hadroniza- 
tion model [6,7], light- and heavy-quark initiated jets show significant differences regarding particle 
multiplicities as a consequence of soft gluon suppression inside the dead cone. Such differences 
could be exploited by using auxiliary criteria complementing ^-tagging procedures to be applied to 
jets with very large transverse momentum, as advocated in this study. 
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2. Kinematics and variables 

As known from jet calculus for light quarks, the evolution time parameter determining the 
structure of the parton branching of the primary gluon is given by (for a review see [4] and refer- 
ences therein) 

y = M7rY k x =zQ>Q 0l Q = E®>Q Q , (2.1) 



where k± is the transverse momentum of the gluon emitted off the light quark, Q is the virtuality 
of the jet (or jet hardness), E the energy of the leading parton, Qo/E < < O is the emission 
angle of the gluon (0 -C 1), 0o the total half opening angle of the jet being fixed by experimental 
requirements, and <2o is the collinear cut-off parameter. Let us define in this context the variable 
Y as y = Y + lnz, Y = In ( -g- ) . The appearance of this scale is a consequence of angular ordering 
(AO) of successive parton branchings in QCD cascades [4, 7]. An important difference in the 
structure of light (£ = q = u,d,s) versus heavy quark (h = Q = c,b) jets stems from the dynamical 
restriction on the phase space of primary gluon radiation in the heavy quark case, where the gluon 
radiation off an energetic quark Q with mass m and energy E S> m is suppressed inside the forward 
cone with an opening angle & m = m/E, the above-mentioned dead cone phenomenon [3]. 

The corresponding evolution time parameter for a jet initiated by a heavy quark with energy 
E and mass m appears in a natural way and reads y = In f ^ J , KJ_ = k\ +z 2 m 2 [3], which for 
collinear emissions 0^1 can also be rewritten in the form 

K± = zQ, Q = E(© 2 + © 2 m y\ (2.2) 

with > 0„, (see Fig. |l|). An additional comment is in order concerning the AO for gluons emitted 
off the heavy quark. In ( |2.2[ ), is the emission angle of the primary gluon g being emitted off 
the heavy quark. Now let 0' be the emission angle of a second gluon g' relative to the primary 
gluon with energy a/ <C CO and 0" the emission angle relative to the heavy quark; in this case the 
incoherence condition /2 < (0 2 + 2 ) together with 0" > m (the emission angle of the second 



gluon should still be larger than the dead cone) naturally leads ( (2.2| ) to become the proper evolution 
parameter for the gluon subjet (for more details see [3]). For m = 0, the standard AO (0' < 0) is 
recovered. Therefore, for a massless quark, the virtuality of the jet simply reduces to Q = E& as 
given above. The same quantity k:^ determines the scale of the running coupling a s in the gluon 
emission off the heavy quark. It can be related to the anomalous dimension of the process by 

where nj is the number of active flavours and N c the number of colours. The variation of the 
effective coupling a s . as nf — > n * + 1 over the heavy quarks threshold has been suggested by next- 
to-leading (NLO) calculations in the MS scheme [6] and is sub-leading in this frame. In this 
context j3o(«/) will be evaluated at the total number of quarks we consider in our application. The 
four scales of the process are related as follows, 

<2>m>£o~A GC Z), 
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where <2o ~ ^qcd corresponds to the limiting spectrum approximation [4]. Finally, the dead cone 
phenomenon imposes the following bounds of integration to the perturbative regime 

g<Z<l-g, m 2 <Q 2 < E 2 (&l + & 2 m ), (2.4) 

which now account for the phase-space of the heavy quark jet. The last inequality states that the 
minimal transverse momentum of the jet Q = E& m = m is given by the mass of the heavy quark, 
which enters the game as the natural cut-off parameter of the perturbative approach. 

3. Definitions and notation 

The multiplicity distribution is defined by the formula 

p n = -p—=*-, £p„ = i (3.1) 

l*n=0 ®n Ginel n= o 

where o n denotes the cross section of an ^-particle yield process, Oi ne i is the inelastic cross-section, 
and the sum runs over all possible values of n. 

It is often more convenient to represent multiplicity distributions by their moments. All such 
sets can be obtained from the generating functional Z(y, u) [4] defined by 

Z(y,«) = £P n (y)(l + «r 

n=0 

at the energy scale y. For fixed y, we can drop this variable from the azimuthally averaged generat- 
ing functional Z(u). The average multiplicity is defined by the formula, 

1 d"Z(u) 



(n)=N='£Pnn, 



n _ {) n\ du n 



(3.2) 
w=-l 

We will compute the average multiplicity ( fO| ) of partons in jets to be denoted hereafter as Na, with 
A = Q,q,g, corresponding to a heavy, light quark or gluon initiated jet respectively. 

Once arrived at his point, let us make an important distinction between two different particle 
sources populating heavy-quark initiated jets. On the one hand, parton cascade from gluon emission 
yields the QCD component of the total jet multiplicity (the main object of our present study), 
excluding weak decay products of the leading primary quark at the final stage of hadronization. 
On the other hand, the latter products coming from the leading flavoured hadron should be taken 
into account in the measured multiplicities of jets. We shall denote the average charged hadron 
multiplicity from the latter source as N^ c . Hence the total charged average multiplicity, N A 0,al , 
reads 

N totai = N ch + N dc . A = q,Q. (3.3) 

As a consequence of the LPHD, N^ 1 = J^ ch x Na [6,7], where the free parameter <%f ch normalizes 
the average multiplicity of partons to the average multiplicity of charged hadrons. For charm and 
bottom quarks, we will respectively set the values N*° = 2.60 ± 0. 15 and N$ c = 5.55 ± 0.09 [3, 8], 
while in light quark jets one expects N? c = 1.2 ±0.1 [9]. In this work, we advocate the use of such 
a difference between average jet multiplicities as a signature to distinguish a posteriori heavy from 
light quark jets, particularly in ^-tagging techniques applied to the analysis of many interesting 
decay channels. 
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4. QCD evolution equations 



The splitting functions 



P{z,a s ) = a s P^(z) + a 2 P {l \z) + ... 



(4.1) 



where P^\z) and P^ x '{z) are respectively the LO and NLO splitting functions, can be associated 
to each vertex of the process in the partonic shower. P(z, oc s ) determines the decay probability of a 
parent parton (quark, anti-quark, gluon) into two offspring partons of energy fractions z and 1 — z. 
In this work, we are rather concerned with calculations which only involve the LO P^ (z) splitting 
functions in the evolution equations [4]. 
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Figure 1: Parton splitting in the process Q 
shown. 



Qg: a dead cone with opening angle ©,„ is schematically 



Let us start by considering the the splitting process, Q — > Qg, Q being a heavy quark and g the 
emitted gluon which is displayed in Fig.|]; the corresponding LO splitting function reads [10, 11] 



Mi 



C F 



Qg W ~ N 



z 2 



21 



z(l -z)m 
kj_ + z 2 m 2 



*-> 



<('--) 



(4.2) 



where k± ~ min(zE®,(l — z)E&) is the transverse momentum of the soft gluon being emitted 
off the heavy quark. The previous formula ( |4.2| ) has the following physical interpretation, for 
kj_ <^i z 2 m 2 , the corresponding limit reads PQ g (z) -^ jfi-z and that is why, at leading logarithmic 
approximation, the forward emission of soft and collinear gluons off the heavy quark becomes 
suppressed once <C ® m , while the emission of hard and collinear gluons dominates in this region. 
For the massless process g — > gg, we adopt the standard three gluon vertex kernel [4, 13] 



P^(z) 



and finally for g — > QQ, we take [10, 1 1] 



(l- z )[2-z(l-z)] ; 



p(°), 

gQ 



)- — 



l-2z(l-z) + 



2z(l — z)m 

k^+m 2 



21 



(4.3) 



(4.4) 



which needs to be resummed together with the three gluon vertex contribution. However, as a first 
approach to this problem, we neglect the production of heavy quark pairs inside gluon and quark 
jets, making use of [4, 13] 

1 



M 



(o), 



W)=Oz)U=o 



AN, 



[l-2z(l-z)]. 



(4.5) 
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Including mass effects in the evolution equations also requires the replacement of the massless 
quark propagator l/kj_ by the massive quark propagator l/(&2 + z 2 m ) [11], such that the phase 
space for soft and collinear gluon emissions off the heavy quark can be written in the form [11] 

(4.6) 



d 2 °Q8^fo ,/ k l^ dzP^(z). 



kj_ + z 2 m 



,(0), 



where PqI(z) is given by (4/2). Working out the structure of (42) and setting kx ~ zE® one has 

1 2 . z ®l 



p (o) () _Cf_ 



& 2 + &l 



1 + 1 + 



2 2 + 02 



(4.7) 



The system of QCD evolution equations for the heavy quark initiated jet is found to read [12] 

^^=e 1 (f,L m )j\y^(y)N g (y)-e2(Y,L m )^(Y)N g (Y), (4.8) 

dN l'^ ev 

~J-= L dy^(y)N g (y)-A(Y,L m )yi(Y)N g (f), (4.9) 

U I J Y,„ 



where 

A(Y,L m )=a(n f ) 
with 



2 + 



2N C 



N C J 



-Y+L„, 



l+^fl-2^ 



a(tif) 



AN r 



>M-m 



N c 



Y =1 

1 in — L ^}i 



N c 



-2Y+2L,,, 



In 



m 



and 



Si(Y,L m ) = l-e 



-2Y+2L m 



e 2 (Y,L n 



3 3 
e 

4 2 



-Y+L,„ 



Qo 

-2Y+2L m 



(4.10) 



(4.11) 



(4.12) 



There are the following kinds of power suppressed corrections to the heavy quark multiplicity: the 



leading integral term of ( |4.8| ) is 0(%) suppressed, while subleading MLLA corrections appear in 

2 

the standard form G^^fa^) like in the massless case, finally ^(2^) and 0(%t/c^), which are 
new in this context. Massless results can be recovered after setting m/Q — > Qo/Q in ( |4.10| ) and 
( ft. 12 ). For massive particles however, these terms are somewhat larger and can not be neglected in 
our approach unless they are evaluated for much higher energies than at present colliders. On top 
of that, the corresponding massless equations in the high energy limit are obtained from ( |4~8| ) and 
|) simply by setting y — > y, Y — > Y, Y^ — > Y, Y m — >■ 0, £\ — Y 1, 



3 3 _ Y 
S2^S 2 = ---e +• 



-2Y 



)> A 



and are written in the standard form [13] 

N c dN, rY 



a(tif) 



2 + 



nj_ 

2N r 



1-2^ 

N c 



e- T + ff{e- 2r ) 



C F dY 



J o dy^N g (y) - h^N g (Y), -£ = J dy^N g (y) -A^N g (Y), (4.13) 



with the initial condition N giq {Y = 0) = 1 at threshold. Notice that ( ft.8[ ) and (4.9) are valid only 
for m S> Qo and therefore m — > does not reproduce the correct limit, which has to be smooth as 
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given by the massless equations (4.13). Since heavy quarks are less sensitive to recoil effects, the 
subtraction terms <=c e +Lm and oc e - 2¥ + 2L m m £ 2 (^£ m ) diminish the role of energy conservation 
as compared to massless quark initiated jets. As a consistency check, upon integration over Y of 



the DLA term in Eq.( |4.8[ ), the phase space structure of the radiated quanta reads 



iV G (lnG)«l + 



C F 



0,1 



& 2 d& 2 



N c Jo (® 2 + & 2 m ) 2 J m /Q z 



\—m 



iQdz 



[$N g ] (InzG). 



(4.14) 



Notice that the lower bound over @ 2 in ( [4. 14[ ) (Y in (p~8|>) can be taken down to "0" (Y m = L m in 



)) because the heavy quark mass plays the role of collinear cut-off parameter. 



5. Phenomenological consequences 



Working out the structure of (4.8) and ( |4.9p , we obtain the rough difference between the light 
and heavy quark jet multiplicities, which yields, 



N n -N c 




N 



q, 




(5.1) 



It can be seen that (5.1) is exponentially increasing because it is dominated by the leading DLA 



energy dependence of N q . According to (5.1), the gap arising from the dead cone effect should be 
bigger for the b than for the c quark at the primary state bremsstrahlung radiation off the heavy 
quark jet. The approximated solution of the evolution equations leads to the rough behaviour of 
N q — Nq in (5A), which is not exact in its present form. In Fig.||| (left), we display the numerical 



solution of the evolution equations ( |4.13| ) for N q and (|4.8| ) for Nq corresponding to the heavy quark 
mass intervals m c = 1 — 1.5 GeV, m& = 3 — 5 GeV. Let us remark that the gap arising between 
the light quark jet multiplicity and the heavy quark jet multiplicity follows the trends given by 
(5.1) asymptotically with E — y °°. Finally, as expected for massless quarks L m = 0, the difference 



N q —Nq vanishes. In this study we advocate the role of mean multiplicities of jets as a potentially 





Figure 2: Massless and massive quark jet average multiplicity Nq for the primary state radiation (left) and 



total 



Nq" c " including heavy quark flavour decays (right) as a function of the jet hardness Q. Bands indicate m c 
and mi, in the [1, 1.5] and [3,5] GeV intervals respectively. 
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useful signature for ^-tagging and new associated physics [14] when combined with other selection 
criteria. In Fig. § (right), we plot as function of the jet hardness Q l , the total average jet multiplicity 
(|Q|), which accounts for the primary state radiation off the heavy quark together with the decay 
products from the final-state flavoured hadrons, which were introduced in section ||. For these 



decays constants N^ c = 2.60 ±0.15 and N^ c = 5.55 ±0.09 are independent of the hard process 
inside the cascade, such that N^ c can be added in the whole energy range. For instance, such 



predictions, we set J^ ch = 0.6 in (3.3), and <2o ~ ^qcd = 230 MeV [4]. Moreover, the flavour 

values were obtained by the OPAL collaboration at the Z° peak of the e + e~ annihilation. In this 
experiment, D* mesons were properly reconstructed in order to provide samples of events with 
varying c and b purity, such that it became possible to measure light and heavy quark charged 
hadron multiplicities separately [8]. As compared to the average multiplicities of the primary state 
radiation displayed in Fig. [| (left), after accounting for N^ c , the b quark jet multiplicity becomes 
slightly higher than the c quark jet multiplicity, although both remain suppressed because of the 
dead cone effect. 

6. Conclusions 

Thus, our present work focusing on the differences of the average charged hadron multiplicity 
between jets initiated by gluons, light or heavy quarks could indeed represent a helpful auxiliary 
criterion to tag such heavy flavours from background for jet hardness Q > 40 GeV. Notice that we 
are suggesting as a potential signature the a posteriori comparison between average jet multiplic- 
ities corresponding to different samples of events where other criteria to discriminate heavy from 
light quark initiated jets were first applied. Fig.|| (right) plainly demonstrate that the separation 
between light quark jets and heavy quark jets is allowed above a few tens of GeV with the fore- 
seen errors of the experimentally measured average multiplicities of jets. The difference between 
light quark jet multiplicities and heavy quark jet multiplicities N q — Nq in one jet is exponentially 
increasing because of suppression of forward gluons in the angular region around the heavy quark 
direction. This result is not drastically affected after accounting for heavy flavour decays multi- 
plicities, such that it can still be used as an important signature for ^-tagging and the search of 
new physics in a jet together with other selection criteria. Notice that the measurement of such 
observables require the previous reconstruction of jets at hadron colliders [15]. 
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